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THE  DELAY-DOPPLER  SPECTRAL  ANALYSIS  OF 
NONSTATIONARY  TIME  VARIANT  LINEAR  CHANNELS 
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I.  INTRODUCTION 
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The  channel  characterization  and  measurement  problem  for 
physical  time  variant  channels  is  complicated  by  nonstationary 
statistical  behavior.  This  complication  is,  in  part,  analogous 
to  that  perceived  in  the  analysis  of  speech  signals,  particularly 
with  a view  to  the  characterization  of  such  signals  in  identi- 
fiable classes.  The  earliest  and  still  popular  approach  to 
speech  analysis  involves  a short-term  running  spectral  analysis 
in  order  to  follow  the  nonstationary  changes  and  to  provide 
identifiable  characteristics  (e.g.,  formants).  A similar  approach 
is  necessary  in  the  measurement  of  nonstationary  channels  and  is 
developed  in  some  detail  in  the  discussion  below.  In  brief,  a 
short-term  running  two-dimensional  delay-Doppler  spectral  analy- 
sis is  studied  as  a means  of  characterizing  the  channel  on  a 
short-term  basis  and  allowing  observance  of  nonstationary  behavior. 
The  analysis  draws  upon  prior  theory  developed  in  [1]  and  [2]  f 

on  the  characterizations  and  measurement  of  randomly  time  variant 

\ 

linear  channels. 
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II.  THE  SPECTRAL  REPRESENTATION  OF 
TIME  VARIANT  CHANNELS 

A variety  of  system  functions  exist  for  characterizing  the 
input-output  relations  of  time  variant  linear  channels.  We  use 
the  categorization  and  definition  of  system  functions  developed 
in  [1],  where  system  functions  are  arranged  in  time-frequency 
dual  pairs.  To  simplify  the  discussion  we  shall  deal  primarily 
with  two  system  functions  T(f,t)  and  U(5,v)  called,  respectively, 
the  time  variant  transfer  function  and  the  delay-Doppler  spread 
function.*  If  z(t)  is  the  complex  envelope  of  the  channel  input 


and  w(t)  is  the  ccmplex  envelope  of  the  channel  output, 
T(f,t)  and  U(£,v)  provide  the  input-output  relations, 

then 

w(t)  = I 
« 

' Z(f)T(f,t)eJ2TTftdf 

(1) 

w(t)  = 

t 

|J  z(t-§)eJ2TTvt  U( 5 , v)d|dv 

(2) 

where  Z(f)  is  the  spectrum  of  z(t). 

The  system  functions  T(f,t)  and  U(|,v)  are  double  Fourier 
transform  pairs 


. J2TTf5-j2TTVt 

(3) 

U(5.v)  = JJ 

T(f ,t)e  e dfdt 

■ -j2n?f  J2 TTvt 

(4) 

T(f,t)  * JJ 

U(5,v)e  e d«;dv 

*The  entire  paper  can  be  phrased  in  terms  of  the  dual  system 
functions  M(t,f)  and  V(v,{)  defined  in  [1], 
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An  examination  of  Eq.  (2)  reveals  that  the  delay-Doppler 
spread  function  provides  a phenomenological  model  of  the  channel 
as  a continuum  of  differential  "scatterers"  subjecting  the  trans- 
mitted signal  to  a complex  gain  U(?,v)d?dv  for  the  scatterers 
providing  delays  (?)  and  Doppler  shifts  (v)  in  the  region 
(?,£+d?)  x (v,v+dv).  The  function  U(?,v)  provides  the  formal 
spectral  representations  of  the  time  variant  channel  that  we 
shall  use  in  this  paper.  It  will  be  used  in  the  same  way  that 
the  Fourier  transform  of  a time  function  is  used  to  provide  a 
formal  spectral  representation  of  a time  series.  It  should  be 
noted  that  we  do  not  propose  to  present  a rigorous  mathematical 
exposition  of  the  Fourier  transform  of  random  processes  but 
rather  an  engineering  presentation  following  the  approach 
developed  by  the  author  in  [3],  Section  III.  Thus  we  shall  side- 
step questions  concerning  the  existence  of  integrals  and  appeal 
to  physical  interpretation  for  the  justification  of  results. 

The  spectral  occupancy  of  a time  function  is  defined  as  the 
set  of  frequencies  over  which  the  spectrum  of  the  function  is 
"significantly"  different  from  zero.  The  bandwidth  of  a time 
function  is  generally  defined  as  the  highest  frequency  of  its 
spectral  occupancy.  However,  in  the  case  of  narrow  band  pro- 
cesses, it  is  more  meaningful  to  define  a "center"  frequency  for 
the  spectral  occupancy  and  then  define  bandwidth  as  the  "width"  of 
the  spectral  occupancy,  i.e.,  the  difference  between  the  highest 
and  lowest  frequencies  in  the  spectral  occupancy. 

For  a stationary  random  process  z(t)  with  power  spectrum 

P (f)  and  autocorrelation  function  R (t),  the  spectrum  of  z(t), 
z z 

defined  as 

r -J2nft 

Z(f)  = J z(t)e  dt  (5) 

has  the  autocorrelation  function  (see  [3],  Section  III), 

Z*( f )Z(4)  - P.(f)«(t-4)  (6) 

z 
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where  6(.)  is  the  unit  impulse  and 

. -j2TTfT 

P (f)  = R ( T)e  dT  (7) 

Z j z 

z*( t)z(  t+r)  = R_(t)  (8) 

z 

The  spectral  occupancy  of  z(t)  will  clearly  not  include 

any  frequencies  for  which  P (f)  is  identically  zero.  Where 

z 

P (f)  is  not  identically  zero  some  criterion  of  "significantly 
z 

different  from  zero"  must  be  established.  When  z(t)  is  a non- 
stationary random  process  the  power  spectrum  does  not  have  any 
meaning  since  the  latter  is  defined  only  for  a stationary  random 
process.  However  the  concepts  of  spectral  occupancy,  bandwidth, 
and  center  frequency  can  still  have  meaning  for  some  nonstationary 
processes. 

Turning  now  to  time-varying  channels  we  define  by  analogy  > 

the  delav-DoPbler  occupancy  pattern  or  just  delay-Doppler 

occupancy  of  the  channel  as  the  set  of  delays  and  Doppler  shifts 

over  which  the  two-dimensional  spectrum  U(£,v)  is  "significantly" 

different  from  zero.  Examples  of  possible  occupancy  patterns 

are  shown  in  Fig.  1 (repeated  from  Fig.  1 of  [2]).  The  parameters 

B and  L denoted  on  these  patterns  are  analogous  to 
ma  x max 

bandwidth  for  a narrow  band  process.  B , the  Doppler  spread, 

[Da  X 

is  the  difference  between  the  maximum  and  minimum  Doppler  shifts 

and  L , the  delay  spread,  is  the  difference  between  the  maxi- 
max ’ 

mum  and  minimum  delays  in  the  occupancy  pattern. 

The  analog  of  the  stationary  process  z(t)  is  the  WSSUS 
(wide  sense  stationary  uncorrelated)  channel  defined  in  [1]  for 
which 

T*(f,t)T(f+n,  t+T)  = R(  0,  t)  (9) 

1 i : i . . ,!  . « . . i • ■ . v I 
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Note  that  the  two  dimensional  correlation  function  of  the  time 
variant  transfer  function  depends  only  on  time  and  frequency 
differences.  It  is  readily  shown  (see  [l]  that  for  the  WSSUS 
channel 

u*( 5, v)u( s+n,  v+u)  = s(5,v)  6(n-5)6(p-v)  (10) 

where  the  two  dimensional  power  density  spectrum  S(§,v)  is  called 
the  scattering  function  [1].  The  delay-Doppler  occupancy  of  the 
channel  will  clearly  not  include  any  (C,v)  values  for  which 
S(5,v)  vanishes.  Where  S(«;,v)  is  not  identically  zero  some 
criterion  of  "significantly  different  from  zero"  must  be  estab- 
lished. When  the  channel  is  nonstationary  the  scattering  function 
has  no  meaning.  However  the  concepts  of  delay-Doppler  occupancy, 
Doppler  spread,  and  delay  spread  can  still  have  significance  for 
some  nonstationary  channels. 

While  the  spectral  representation  of  a channel  is  a clear 
generalization  to  two  dimensions  of  the  spectral  representation 
of  a time  function,  this  is  not  true  for  the  spectral  measurement  of 
a time  series.  We  note  that  in  the  channel  measurement  case 
one  must  transmit  probing  waveforms  through  the  channel  to  gener- 
ate signals  related  to  the  time  variant  transfer  function.  In 
addition,  under  certain  conditions  the  channel  system  functions 
are  nonmeasurable  even  in  the  absence  of  additive  noise.  Thus 
the  author  has  shown  [2]  that  a time  variant  channel  is  measurable 
if  and  only  if  the  area  of  the  delay-Doppler  occupancy  pattern 
is  less  than  unity.  Finally,  the  measurement  problem  in  the 
presence  of  noise  [2]  is  considerably  more  difficult  to  deal 
with  in  the  channel  measurement  case  even  apart  from  the  Increase 
in  dimensionality  of  the  problem. 
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III.  SHORT-TERM  SPECTRAL  ANALYSIS  OF 
TIME  VARIANT  CHANNELS 


The  integrals  in  (1)  and  (3)  involve  integration  over  all 
time  and  frequency  and  are  thus  convenient  mathematical 
fictions.  Although  only  finite  records  of  a process  in  time 
and  frequency  can  ever  be  available  and  although  statistical 
stationarity  is  the  exception  rather  than  the  rule,  this  has  not 
prevented  spectral  analysis  from  being  an  extremely  useful  tool 
to  the  communication  engineer  in  the  characterization  and  pro- 
cessing of  communication  signals  and  interference.  The  type  of 
spectral  analysis  that  is  employed  in  practice  is  based  upon 
finite  intervals  of  time  and  sometimes  short-term  running 
intervals  of  time.  A case  in  point  is  the  analysis  of  speech 
signals  with  a Sonogram  [3]  which  presents  a short-term  intensity 
spectrum  as  a function  of  time.  The  advantage  of  such  a running 
analysis  is  that  "quasi-stationary"  spectral  characteristics 
can  be  observed.  The  purpose  of  this  section  is  to  define  such 
a short  term  delay-Doppler  spectral  analysis  for  time  variant 
channels  and  discuss  some  of  its  properties. 


3.1  The  Segmented  Channel 


Just  as  the  basis  for  the  short-term  spectral  analysis  of 
an  analog  waveform  is  the  Fourier  analysis  of  a finite  segment 
of  the  waveform,  the  basis  for  the  short-term  spectral  analysis 
of  a time  variant  channel  is  the  Fourier  analysis  of  a "finite 
segment"  of  the  channel.  Here  "finite  segment"  is  interpreted 
as  a hypothetical  channel  whose  time  variant  transfer  function  is 
identical  to  the  transfer  function  T(f,t)  of  the  actual  channel 
in  a rectangular  region  of  t-f  space  and  is  zero  elsewhere.  Thus 
we  define  a truncated  time  variant  transfer  function 


where 


T(f,t)  ■ Rect(f/W)T(f , t)Rect(tA) 


r 

1 


Rect(x) 


* 
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Nl  < i 
1*1  i | 


(11) 


(12) 
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One  may  obtain  a physical  representation  of  a channel 

A 

with  time  variant  transfer  function  T(f,t)  in  terms  of  the 
channel  with  time  variant  transfer  function  T(f,t)  by  preceding 
the  latter  with  a filter  having  transfer  function  Rect(f/W) 
(i.e.^an  ideal  band  pass  filter  of  bandwidth  W Hz)  and  following 
it  with  a time  jate  Rect(t/T)  of  duration  T seconds.  Using  the 
Fourier  transform  property  (3),  (4)  it  is  readily  shown  that  the 
spectral  representation  of  the  segmented  channel  is  related  to 
that  of  the  original  channel  by 


where 


U(?,v)  = JJ  TW  sine  [W(  rr?)  ] T sine  [T(u-v)  ]U(  n,p)  di-idp 

(13) 


sine  x 


sin  ttx 

TTX 


(14) 


The  spectrum  of  the  segmented  channel  is  obtained  by  a two 
dimensional  convolution  of  the  spectrum  of  the  original  channel 
with  the  spectrum  TW  sinc(W§)sinc(Tv) . It  is  important  to  note 
that  this  latter  spectrum  is  non-vanishing  in  any  finite  non- 
zero interval  of  5 or  v and  decreases  very  slowly  from  its 
maximum  value  of  unity.  Thus  although  the  delay-Doppler  occu- 
pancy pattern  of  the  original  channel  may  be  identically  zero 
outside  a given  region  of  the  ^,v  plane.  The  occupancy  pattern 

of  the  segmented  channel  cannot  be  identically  zero  over  any 

* 

nondegenerate  region.  Of  course,  the  occupancy  pattern  of  the 
segmented  channel  will  become  small  sufficiently  far  from  the 
occupancy  pattern  of  the  original  channel  and  will  eventually 
become  sufficiently  small  to  be  regarded  as  negligible.  How- 
ever the  occupancy  pattern  will  be  a spread-out  version  of  the 
original,  the  amount  of  spreading  decreasing  as  T and  W increase. 

A discrete  spectral  representation  of  the  segmented  channel 
is  provided  by  the  use  of  sampling  theory  applied  to  (5(5,^). 
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Thus  since  the.  transform  of  U(g,n)»  T(f,t),  is  "band-limited" 
in  both  the  variables  t and  f we  may  reconstruct  U(§,h)  from 
the  samples 

Umn  S TO  u(w  » t)  ^15' 

A 

The  impulse  sampled  version  of  U(|,r|)  is  given  by 

- 12  °mn6(5-  §)»(v-  |)  (16! 

m n 

where,  from  (15)  and  (13), 

Umn  = If  sinc  LW(§  “ w)j  sinc  [_T(v-?)jU(5,v)d5dn 


The  quantity  Umn  is  readily  shown  to  be  equal  also  to  the  generic 
coefficient  in  the  Fourier  expansion  of  T(f,t)  (and  T(f,t)) 


T/2  »W/2 

TO  I J T(  f , t)  e 
-T/2  -W/2 


J2TTf  g-  J2Tr|t 


df  dt 


The  segmented  channel  is  "recovered"  from  the  channel  which 

has  the  impulse  lattice  (16)  for  its  delay-Doppler  spread  function 

by  the  use  of  band  limiting  before  and  time  gating  after  the 

latter  channel  as  shown  in  Fig.  2.  This  is  the  discrete  delay- 

Doppler  representation  discussed  in  [1]  and  [2],  The  input -output 

relation  of  the  segmented  channel  is  uniquely  determined  by  the 

coefficients  fU  }.  Thus  if  the  segmented  channel  input  z(t)  is 
mn 

limited  to  the  band  of  frequencies  -W/2  < f < W/2,  the  output  w(t) 
is  given  by 


j 2 TT  “ ■ t 

w(t)  = Rect( t/T)  ZE  U z(t  - £j)e  (19) 

mn  mn  W 

If  the  input  is  not  bandlimited, then  z(*)  in  (19)  should  be 
replaced  by  only  those  frequency  components  of  z(.)  which  lie 
within -W/2  < f < W/2. 
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The  number  of  U coefficients  which  are  significantly 
ran 

different  from  zero  depends  upon  the  occupancy  pattern  of 
U(^,v)  and  the  values  of  T and  W as  illustrated  in  Fig.  3. 

Also  shown  in  Fig.  3 is  the  occupancy  pattern  of  the  original 
unsegmented  channel  in  addition  to  that  of  the  segmented 
channel,  showing  the  spreading  of  the  latter  relative  to  the 
former.  The  spreading  of  the  occupancy  pattern  after  segmenta- 
tion is  unfortunate  because  it  increases  the  number  of  spectral 
coefficients  that  need  to  be  measured  in  order  to  cnaracterize 
the  segmented  channel.  In  Section  3.3  we  consider  the  use  of 
window  functions  as  a means  of  reducing  the  spreading  of  the 
occupancy  pattern  of  segmented  channels. 

3.2  Representation  of  Channel  in  Terms  of  Segmented  Channels 

The  above  discussion  was  confined  to  a segmented  version 
of  the  actual  channel  which  accepted  frequencies  in  the  band 
-W/2  < f < W/2  and  had  output  only  in  the  time  interval 
-T/2  < t < T/2.  To  provide  spectral  estimates  over  wider  fre- 
quency and  time  intervals  one  may  either  increase  T and  W or 
else  keep  T and  W fixed  and  use  more  segmented  channels  to 
characterize  the  additional  time  and  frequency  intervals.  In 
order  to  observe  nonstationary  characteristics  in  time  and  fre- 
quency one  should  keep  T and  W as  small  as  possible  and  build 
up  the  actual  channel  with  segmented  channels. 

Thus  we  consider  a channel  with  time-variant  transfer 
function  given  by  the  series 

A P Q 

T(f,t)  - E E T (f , t)  (20) 

p=_p  q=_Q  P<3 

where  T (f,t)  is  the  time  variant  transfer  function  of  a com- 

pq 

ponent  segmented  channel 

Tpq(f,t)  = Rect(|  - p)T(f,t)Rect(|  - q)  (21) 
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in  which  the  input  frequency  band  is  confined  to  pW  - W/2  < f 
< pW  + W/2  and  the  output  time  to  qT  - T/2  < t < qT  + T/2. 


Note  that  since 


E Rect(-  - p)  = Rect[(2P  + l)W_] 


p=-P 


Z Rect (t  “ q)  Rect[(2Q  + 1 )tJ 
q=  -Q 


(22) 


the  summation  (19)  is  identical  to 


T(fft)  = RectQ- 


2P  + 1)WJ 


T(f 


, t)Rect  ■ 


(2Q  + 1 ) T 


J 


(23) 


which  is  a representation  of  the  channel  for  input  frequencies 
-(P  + i/2)W  < f < (P  + 1/2)W  and  output  time  instants 
-(Q  + 1/2)T  < t < (Q  + 1/2)T. 

By  analogy  with  Eq.  (13),  the  delay-Doppler  spread  function 
of  the  component  segmented  channel  (20)  is  readily  found  to  be 


U <5,v)  = JJ  e-J2TTpW(T1-5)eJ2TTqT(u"v) 

P<1  , (24) 

TW  sine  [W(n-5)  j sxnc  [T(u-v)  ]U(  r|,|j)dr|du 

The  impulse  sampled  version  of  U^t^v)  analogous  to  (16)  is 


0 (5,v)  * i ruM  «({-§)  «(?  -f) 

m n 


(25) 
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where  the  spectral  samples  are 


mn  TW  pq 


Vw  » T ) 


= JJ  e_j2rrpW5ej2TTqTv  sine  [w(?  - g)j  sine  [t(v  - |)j 

U(5,v)d§dv  (26) 

As  in  (18)  the  spectral  samples  are  expressible  as  Fourier 


coefficients 


nrT  .qT+T/2  pW  + W/2  J2irf  g -J2iT§t 

"S-*;  i T(f-t,e  • a£at 


qT  - T/2  pW  - W/2 


The  set  of  spectral  samples  for  fixed  values  of  p 

and  q provide  a "snapshot"  spectral  analysis  of  the  channel. 

As  p and  q are  changed  one  obtains  a discrete  "running"  spectral 
analysis  in  the  (output)  time  and  (input)  frequency  variables. 

3 . 3 Window  Functions 

In  practice  only  a finite  number  of  spectral  samples  can 
be  used  to  characterize  a component  segmented  channel.  As  dis- 
cussed at  the  end  of  Section  3.1  and  illustrated  in  Fig.  3, 
the  number  of  spectral  samples  needed  is  determined  by  the  area 
of  the  occupancy  pattern  of  the  component  segmented  channel. 

This  area  in  turn  is  dependent  upon  some  arbitrary  criterion  as 
to  what  values  of  Upq(5,v)  are  "significantly"  different  from 
zero.  In  the  discussion  that  follows  we  assume  that  a threshold 

criterion  has  been  established  which  allows  us  to  regard  U (t,v) 

pq  »» 

as  negligible  when  it  is  less  than  some  constant  c.  Thus  if 

A 

Sp<^  denotes  the  set  of  points  in  the  occupancy  pattern  of  the 
segmented  channel, 
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(?,v)  eSpq  s=^|Upq(?,v)|  > c Max|Upq(5fv)|  (20) 

The  occupancy  pattern  of  the  segmented  channel  is  a somewhat 
spread-out  version  of  the  actual  channel  occupancy  pattern.  Thus 
in  the  case  of  a channel  consisting  of  a single  delay  and  Doppler 
shift,  i.e., 

U(S,v)  = - 5j)Mv  - v:)  (29) 

we  note  from  Eq.  (23)  that  the  magnitude  of  the  spectrum  of  the 
segmented  channel  is  given  by 

|Upq(§,v)|  = l^TW  sine  [W(51  - 5)]  sine  [T^  - v)  ] 

(30) 

Since  the  sine  function  drops  to  zero  very  slowly  and  in  an 
oscillatory  fashion,  the  zero-area  occupancy  pattern  of  the  original 
channel  expands  into  a pattern  of  occupied  regions.  For  example, 
if  c in  (27)  is  defined  as  10%,  the  occupancy  pattern  of  the  seg- 
mented channel  takes  the  form  shown  in  solid  lines  in  Fig.  4. 

While  the  total  area  occupied  by  this  pattern  is  approximately 
4.5/IW,  note  that  it  extends  almost  to  +3/W  in  5 and  +3/T  in  v. 

When  c is  decreased  to  3%  the  pattern  broadens  drastically  as  shown 
in  Fig.  5.  The  area  has  increased  to  approximately  27 ,8/TW  and 
the  pattern  extends  to  +8/W  in  5 and  +8/T  in  v. 

In  order  to  observe  the  (time  and  frequency)  non-stationary 
characteristics  of  the  channel  it  is  desirable  to  make  T and  W as 
small  as  possible.  However  we  see  that  the  penalty  is  an  increase 
in  the  area  of  the  segmented  channel's  occupancy  pattern  and  a 
consequent  increase  in  the  number  of  discrete  spectral  samples 
needed  to  characterize  each  segmented  channel.  Note  that  if  TW 
becomes  small  enough  the  segmented  channel  could  become  non- 
measurablel 
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In  order  to  alleviate  this  spreading  problem  one  may  use  an 
expansion  of  the  time  variant  transfer  function  over  the  fre- 
quency band  -(P  + 1/2)W  < f < (p  + 1/2)W  and  time  interval 
-(Q  + 1/2 )T  < t < (Q  + ]/2)T  which  is  different  from  that  shown 
in  Eq.  (19)  in  that  rectangular  time  and  frequency  gates  are  not 
used. 


Consider  the  representation 

P+a  Q+p 

fl(f,t)  = I E 

p=-( P+a)  q=-(Q+B) 


- pW ) G ( t - qT)T(f,t) 


where 


E H(f  - pW)  = 1 


E G(t  - qT)  = 1 (33) 

— ao 

H(f)  = 0 ; \f 1 > aW/2  (34) 

G(t)  =0  ; | t 1 > 0T/2  (35) 

We  call  any  function  which  satisfies  (32)  and  (34)  or  (33)  and  (35) 
a window  function  of  bandwidth  aw  or  time  duration  0T,  respectively. 
Nyquist  (.[4]  has  shown  that  large  classes  of  functions  can  be 
generated  to  satisfy  what  we  call  the  window  property  by  starting 
with  a rectangle  and  adding  odd  symmetric  functions  at  the  edges, 
keeping  the  resulting  function  symmetrical.  Typical  normalized 
(W,T  = 1)  window  functions  arec(*)  the  "raised-cosine"  function 
and  Tr(.)  the  triangle  function,  defined  by 


C(x)  = 


2+2  COS 


Ixl  < 1 


•I  * 1 


Tr(x)  = 


1 - |x|  ; |x|  < l 

0 ; Ixj  > 1 


For  purposes  of  discussion  let  us  use 

H(f)  = C( f/W) 

G(t)  = C( t/T) 


Note  that 


and  that 


a = B = 2 


P+1  Q+l  ,f  «.  ,y  \ 

I I C(^  - p)  = q) 

P=-(P+1)  q=-(Q+l)  K ' 


-( P+1/2 )W  < f < ( P+1/2 )W 
-(Q+l/2 )T  < t < ( Q+l /2 ) T 


A 1 Q+ » If  \ if  \ 

T(f,t)  = I X c(h  - p)  c(£  - q)  T(f,t) 

p=-(P+l)  q=-(Q+l ) XW  ' VT  ' 


is  identical  to  T(f,t)  over  the  input  frequency  range  -(P+1/2)W 
< f < ( P+1 /2 )W  and  the  output  time  interval  -(Q+1/2)T  < t < (Q+1/2)T, 
If  we  define 


Thus 


W£.'>  ' c(s  - p)  c(i  - <0 
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then  T(f,t)  is  represented  as  a summation  of  component  segmented 
channels  with  time  and  frequency  weightings  applied  by  the  window 
functions.  Equation  (41)  may  be  compared  to  Eq.  (20)  which  can 
now  be  viewed  as  the  case  of  rectangular  window  functions.  It  is 
important  to  note  that  the  duration  of  the  windows  in  (41)  are 
2W  and  2T  as  opposed  to  W and  T in  (20).  Since  adjacent  fre- 
quency windows  are  separated  by  W Hz  and  adjacent  time  windows 
by  T sec,  the  windows  overlap  in  (40)  and  not  in  (19).  Also  note 
that  the  summation  in  (40)  extends  one  unit  farther  in  each  direc- 
tion compared  to  (19)  so  that  somewhat  more  component  segmented 
channels  are  required. 

These  additional  complexities  are  offset  by  a reduced  spread- 
ing of  the  occupancy  pattern  of  the  segmented  channel.  Thus 
upon  Fourier  transforming  (41),  we  note  that  the  spectrum  of  the 
segmented  channel  with  raised  cosine  windows  is  given  by 


U (?,v)  = JJ  e-J2T,PW(Tl“5)eJ2T,qT(u“v>TW 


2 * ' 2 U(ri#n)dTidu 
At.  . . \ ^ \ 


(1  - 4W2(t1  - S)2)(l  " 4T2U  “ v)2) 


where 

c(y)  = iLknc_2y  (44) 

1 - 4y^ 

is  the  Fourier  transform  of  the  raised  cosine  window  function  C(x) 

To  see  the  reduced  spreading  of  the  segmented  channel's 
occupancy  pattern  associated  with  raised  cosine  windows  we  con- 
sider the  case  of  a channel  consisting  of  a single  delay  and 
Doppler  shift,  Eq.  (28).  Using  (28)  in  (42)  we  find  that  for 
raised  cosine  weighting 


18 


|Upq( ^,v) \ sine  [2W(g1  - Q]  sine  [2T(v1  - v)  ] 

Max  |Upq(5,v>!  ’ [1  - 4W2(5l  - 5)2][1  - 4T2(Vl  - v)2] 

(45) 

For  c = 10%  the  occupancy  pattern  takes  the  form  shown  by 
dashed  lines  in  Fig.  4.  The  area  of  this  pattern  is  2.3/TW  and 
extends  to  around  +1/W  in  5 and  +1/T  in  v.  The  improvement  of 
raised  cosine  weighting  over  rectangular  weighting  becomes  very 
noticeable  for  c = 3%  as  shown  in  Fig.  5.  Here  the  area  has 
increased  to  only  3.6/TW  and  the  pattern  still  extends  to  only 
around  +1/W  in  5 and  +1/T  in  v. 

For  general  window  functions  it  is  clear  that  the  spectrum 
of  the  segmented  channel  can  be  expressed  as 


Upq<?iV)  = JJ  e-J2"pW(r'-!)eJ2^T(“-'',hln-5)g(U-v)Upq(^u)<S^u 

(46) 

where  h(£),  g(v)  are  the  Fourier  transforms  of  the  window  functions 
H(f)  and  G(t),  respectively. 

Consider  now  the  discrete  representation  of  each  of  the  seg- 
mented channels.  We  note  first  that  the  application  of  general 

window  functions  as  described  above  causes  an  increase  in  the  input 

• • \ 

bandwidth  and  output  time  duration  o*  r which  the  segmented  channel 
is  characterized  to  aW  and  BT,  respectively.  The  sampling  grid 
then  becomes  finer  with  samples  occuring  at  (m/aW,  n/0T)  where  m,n 
are  integers.  For  the  case  of  the  raised  cosine  windows  the  samples 
are  at  (m/2W,  n/2T) . Thus,  in  the  latter  case,  the  discrete  spectral 
samples  needed  to  represent  the  channel  are  given  by 
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I 


0 


u 


pq 

ran 


1 

4TW 


( J (JL  JL\ 
pq\2W  , 2T/ 


■ 4 If 


e - jTTpme  j TTqne  - j 2 TipW§  e j 2 nqTv 


(47) 


- sine  2W(c  - m/2W)  sine  2T(v  - n/2T)  . 

[1  - 4WZ(?  - m/2Wr][l  - 2TZ  ( v - n/2T)^] 

It  should  be  clear  that  the  number  of  discrete  spectral  coefficients 
required  to  characterize  the  segmented  channel  with  raised  cosine 
windows  of  durations  2T,  2W  is  around  four  times  that  required  for 
rectangular  windows  of  durations  T,  W,  at  least  ignoring  the  differ- 
ences of  spreading  of  the  occupancy  patterns  for  the  two  cases. 

One  may  use  raised  cosine  windows  of  smaller  width  than  2T,  2W, 
but  then  more  segmented  channels  would  be  needed  to  characterize 
the  channel  over  large  intervals  of  input  bandwidth  and  output  time. 

For  example,  suppose  the  raised  cosine  windows  were  of  duration  T,  W. 
Roughly  four  times  as  many  segmented  channels  would  be  required  but 
each  segmented  channel  would  require  roughly  the  same  number  of 
spectral  samples  as  rectangular  gates  of  duration  T,  W if  the  oc- 
cupancy patterns  of  the  segmented  channels  had  the  same  area  in  both 
cases.  However,  as  Fig.  5 shows,  the  spreading  caused  by  the  rec- 
tangular gates  can  be  rather  severe  as  compared  to  raised  cosine 
gates. 

Assuming  that  the  window  widths  of  2T  and  2W  are  used  and  thus 
the  spectral  samples  are  given  by  (46),  the  input-output  relation 
for  the  segmented  channel  becomes 


wpq(t)  = ReCt  (2T  " f)  ^ £ Umn  zp(t  ” 2w) 
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(48) 


IV.  SPECTRAL  MEASUREMENT 

In  a previous  paper  [2]  the  channel  measurement  problem  was 
formulated  for  a single  segmented  channel  with  rectangular  windows 
on  the  input  time  and  output  frequency.  Here  we  wish  to  consider 
the  more  practical  running  discrete  spectral  analysis  with  shaped 
windows  outlined  in  the  previous  section.  For  concreteness  we  will 
assume  overlapping  raised  cosine  windows  of  duration  2T  and  2W 
for  input  frequency  and  output  time.  We  first  confine  our  atten- 
tion to  the  case  of  a single  input  frequency  window  but  a large 
number  of  overlapping  output  time  windows.  As  shown  in  Fig.  6, 
the  frequency  window  is  assumed  to  cover  the  band  -W  < f < W Hz, 
Also  it  is  assumed  that  the  input  signal  z(t)  contains  no 
frequency  components  outside  of  this  band  and  has  a non-zero 
spectrum  at  all  frequencies  within  this  band,  except  possibly  for 
f = +W.  From  (48)  we  find  that  the  output  of  a particular  time 
window  located  at  (q  - 1)T<  t < (q  + 1)T  may  be  expressed  as 


« (t>  = Rect  (£  - f)  S t U^n  z(t  - §)  e 


. n . 

JtT  jp  t 


m n 


+ c(|  “ q)n(t) 


(49) 


where  n(t)  is  an  additive  noise,  assumed  white  in  the  present  dis- 
cussion. The  "short-term"  spectral  samples  are  given  by 


uq  = ejTrqn 
ran  — t — 


JJ  e^2rTqTv  c(W5  - m)c(Tv  - n)U(^,v)dcdv 

(50) 

where  c( • ) is  the  transform  of  the  raised  cosine  spectrum  (see  (4-31)) 


24 


The  least-squares  measurement  technique  to  extract  the 
spectral  samples  {u^n}  from  w^(t)  can  be  approached  exactly  as 
in  [2],  However,  here  we  follow  a slightly  different  procedure 
using  a continuous  rather  than  discrete  representation  of  w^(t) 
(see,  for  example,  the  formulation  by  Levin  [5]).  To  simplify 
the  notation  we  assume  the  set  of  possible  discrete  delay- 
Doppler-shif t pairs  has  been  ordered  in  some  convenient  fashion. 
The  k delay-Doppler  shift  pair  (m^,  n^)  is  denoted  by  the 
vector  <1  . The  received  waveform  can  then  be  represented  by  the 


vector  product 

w (t)  = Rect  uTz(t)  + c(£  - q)n(t) 


(51) 


T T 

where  the  vectors  U and  Z (t)  are  qiven  by 

~q 


£ - c»5  . 


d-i 


-2 ' 


V 


(52) 


Z (t)  = [z.  (t),  z,  (t),  ...  z (t)] 
~1  ~2  as 


(53) 


in  which 


* CVV1 


= U<3  „ 
^k  mk'nk 


(54) 

(55) 


I"  jTT^-t 

Z0^t)  = z(t  " 2W)e 


th 


(56) 


The  discrete  spectral  representation  of  the  q channel  is  given 
by  the  vector  U^.  The  least  squares  estimate  of  is  given  by 
(see  [5],  pg.  40) 

U = E^M  (57) 

~q  ~q  ~q 
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where  the  typical  element  of  E and  M are  given  by 


j(q+l)T 

(q-l)T 


z,  (t)z*  ( t ) dt 

Ik  ii 


«!  - 1 


(q+l)T 


(q-l)T 


z,  (t)w  (t)dt 
q 


The  matrix  M is  essentially  the  set  of  matched  filter  outputs, 

~q 

One  may  write  (59)  in  the  alternate  form 


,(q+l  )T 

U = E"iZ*(t)w^(t)dt 

-q  J ~q  ~ q 

(q-l)T 


If  we  define  the  column  matrix 


e ( t)  = E~1Z*(t) 

~q  ~q  ~ 


then 


r(q+i)T 

U = e (t)w  (t)dt 

(q-l)T 


The  matrix  e (t)  is  a column  vector  of  s time  functions  which 

~q 

may  be  computed  once  the  occupancy  pattern  of  the  channel  is 
known.  Equation  (62)  represents  the  estimate  of  each  spectral 
coefficient  as  an  inner  product  of  the  received  signal  and  one 
of  the  component  time  functions  of  e (t).  Alternatively,  Eq.  (57) 

~q 

shows  that  this  estimate  may  also  be  represented  as  a weighted 
sum  of  "matched  filter"  outputs. 
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The  least  square  estimate  (57)  is  the  minimum  variance  linear 
estimator  only  when  the  noise  is  white  and  WSS  (wide  sense 
stationary) . For  more  general  noise  the  minimum  variance  linear 
estimate  is  the  Markov  estimate  Tsee  [5]].  In  the  channel  under 
study  we  consider  only  the  case  of  white  and  WSS  noise  at  the 
receiver  input.  However,  it  should  be  noted  that  the  imposition 
of  the  segmenting  time  gate  makes  the  noise  observed  after  the 
time  gate  nonstationary.  Thus  it  might  appear  as  if  the  more 
general  Markov  estimate  should  be  used.  For  our  special  case  we 
may  obtain  the  Markov  estimate  by  the  simple  artifice  of  dividing 
both  sides  of  (49)  or  (51)  by  the  time  gate  C(|j  - g),  which  removes 
the  raised  cosine  gate  from  both  the  noise  and  observed  data, 
replacing  it  by  the  rectangular  gate  Rect(t/2T  - q/2).  Equation 
(51)  becomes 

w(t)Rect(t/2T  - q/2)  = ^(t)  + n(t) 

(63) 


where  w(t)  is  the  receiver  input  waveform.  Since  the  noise  in 
(63)  is  white  and  WSS  now  one  may  employ  least  squares  to  obtain 
the  minimum  variance  linear  estimate.  Thus  if  this  estimate  is 

A A 

denoted  by  U we  find 


A A 

U 

~q 


F-1N 

~q  ~q 


(64) 


where  the  typical  elements  of  F~*  and  N are 


»(q+ 1 )T  zd  (t)zd.(t) 

F° . = J — ^ dt 


k/. 


(q-l)T  |C(tA-q)r 


r(q+l)T  Zd  (t>W(C) 
N4  “ J C*Tt7T“^"qF 


(65) 


dt 


(66) 


1 


An  examination  of  (65)  and  (66)  reveals  that  certain  diffi- 
culties may  occur  in  applying  (64)  since  the  gating  function  can 


drop  to  zero  and  the  coefficients  and  might  become  infinite. 

IS. 


This  singular  case  has  not  been  examined  yet.  The  covariance 
matrix  of  the  measurement  error  for  the  least  squares  and  Markov 
estimates  are  [see  [5]  and  [2]]. 


c (u  ) 

ov'  _q 


AA  A A *T* 

<u  u*  > - <u  > <u*  > 
~q~q  ~q  ~q 


(67) 


Cov(“q>  * ^<ln(t)12> 


(68) 


where  the  typical  element  of  the  matrix  D is  given  by 

~q 


& - 


k/. 


»(q+l  )T 

j |C(tA  - q)  | zd  (t)z*  (t)  dt 

(q-l)T  U 


(69) 


We  consider  now  the  more  general  case  in  which  there  are 
2P+1  input  raised  cosine  frequency  windows  each  of  width  2W  cps, 
with  the  "center"  window  located  at  -W<f<W  Hz.  From  Eq.  (48)  we 
see  that  the  output  of  a particular  time  window  located  at  (q-l)T 
<f<(q+l)T  may  be  expressed  as 


wqtt)  * 1 wpq(t) 

p=-P 


jTTm  t 


= Rect( t/2T  - q/2)  Eli  U^zo( t - ^_)e"  T 


p m n 


mn  p 


+ C(tA  - q)n(t) 


(70) 
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where  z (t)  is  that  portion  of  the  input  signal  z(t)  located 
P 

within  the  spectral  band  (p-1 )W<f <( p+1 )W.  Equation  (51)  may 

still  be  used  to  represent  w (t)  provided  we  used  the  definitions 

^3 


zT(t)  = [z3p(t)|z'J_p(t)|...|z^1(t)lzJ(t)|z'5'(t)|...|z^t)] 


in  which 


T 

U 

~p,q 


[upq  uPq 

L d * d • 
-1  X2 


upq] 

• • » J 


5p(t>  = C'pdj^-  zpd2't*  ’ 


• zpd  (t)l 


where 


Ik  ■ &V"k] 


uPq  _ L-pq 

Vnv 


m.  jTTT^t  (77) 

zPdk(t)  = zP(t  “ 2w)e 

In  the  definitions  above  it  has  been  assumed,  for  simplicity,  that 
the  occupancy  pattern  is  the  same  for  each  segmented  channel.  If 
desired,  one  may  extend  the  notation  to  a different  occupancy 
pattern  for  each  segmented  channel.  The  discrete  spectral  repre- 
sentation of  the  segmented  channel  with  the  pth  frequency  gate 
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I w 


th 

and  q time  gate  is  the  vector  U . The  set  of  spectral  repre- 

~Pq  *-v 

sentations  for  all  the  frequency  gates  and  the  q time  gate  is 

A 

given  by  U The  least  squares  estimate  of  U , U can  be  expressed 
~q  ~ q 

in  the  general  form  (57)  where  U , E and  M must  be  redefined  as 

~q  ~q  ~q 


the  following  partitioned  matrices 


~q  ,_p  , ~p  ~q,_p.  i~p 


| ~q’i—P , -P  ~q ’ l-P , l-P 


F 

•••  rr  »P- 


~q ’P-1 , P-1  I ?q'P-l,P 


**•  ~q  'P , P—1  ^q’P,P 


- ^P,q^-P,ql 


■•1*5 


• '-P.cr^ 


uT  = roT  i UT  ' ^uT  i i uT  i 

1 X— P 9 q * *,1  — p c * " * * io  q ' • • • l ^ p q J 

The  typical  elements  of  E and  M are 

-qP^  HPq 


qpip2  r(q+i)T 

Sk-i  = J 

(q-l)T 


VT  - - 2s>e  T * 4 < 
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j(q+l)T 

(q-l)T 


zj(t 


-j£n.t 

2 W')e  dt 


Since  the  time  functions  z (t)  and  z (t)  occupy  portions 

P1  p2 

of  the  same  frequency  band  only  when  p^  and  p2  differ  by  unity,  as 
long  as  IW  » 1 and  the  Doppler  spread  of  the  channel  is  much  less 
than  W. 


qpip2 


» 0 for 


lPl  - p2l  > 1 


and  the  matrix  will  be  sparse. 

Markov  estimates  and  error  variances  for  the  spectral  coefficients 
may  be  formed  as  in  the  single-frequency  gate  case. 
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